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Abstract 



We calculate the magnetic moment ('persistent current') in a strongly cor- 
related electron system — a Wigner crystal — in a one-dimensional ballistic 
ring. The flux and temperature dependence of the persistent current in a per- 
fect ring is shown to be essentially the same as for a system of non-interacting 
electrons. In contrast, by incorporating into the ring geometry a tunnel barrier 
that pins the Wigner crystal, the current is suppressed and its temperature 
dependence is drastically changed. The competition between two tempera- 
ture effects — the reduced barrier height for macroscopic tunneling and loss 
of quantum coherence — may result in a sharp peak in the temperature de- 
pendence. The character of the macroscopic quantum tunneling of a Wigner 
crystal ring is dictated by the strength of pinning. At strong pinning the 
tunneling of a rigid Wigner crystal chain is highly inhomogeneous, and the 
persistent current has a well-defined peak at T ~ 0.5 UsjL independent of the 
barrier height (s is the sound velocity of the Wigner crystal, L is the length of 
the ring). In the weak pinning regime, the Wigner crystal tunnels through the 
barrier as a whole and if > Tq the effect of the barrier is to suppress the cur- 
rent amplitude and to shift the crossover temperature from To to T* ~ ^VpTo. 
{Vp is the amplitude of the pinning potential, Tq = Hvp/L, vp ~ h/ma is the 
drift velocity of a Wigner crystal ring with lattice spacing a) . For very weak 
pinning, Vp <^Tq, the influence of the barrier on the persistent current of a 

Wigner crystal ring is negligibly small. 
PACS numbers: 73.40.Gk, 7340.Jn 
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I. INTRODUCTION 



In multiply-connected condensed matter systems, the Aharonov-Bolim effect manifests 
itself as oscillations of both thermodynamic and kinetic quantities with respect to the en- 
closed magnetic flux, $ (for a review see Ref. |l]). From a purely theoretical viewpoint, the 
most evident consequence of these oscillations is the persistent current flow /($) = —cdF/d^ 
(here F is the free energy) induced by the magnetic flux that penetrates through quasi-lD 
conducting ringsi. In normal metal-rings the oscillation period is exactly the quantum flux 
unit $0 = hc/e, where e is the electron charge, and for fixed circumference L, the oscillation 
amplitude can not exceed Iq ~ evp/L. The limiting value of the current, Jq, is determined 
by the diamagnetic electron flow in a Fermi surface orbital and can be evaluated! from a 
quite naive model considering the free electron gas in a perfect ring (i.e. a ring without any 
impurities; see also the review in Ref. 

While the theoretical prediction of such persistent current flow in normal metal-ringsi 
was made some time ago, the corresponding experiments were performed only in the late 
80 's. The measurement of the persistent current was reported in Ref. |^, where the induced 
magnetic moment in a stack of 10^ copper-rings was detected. The period of oscillations 
was found to be $o/2, rather than $o- The apparent puzzle was resolved by noting that the 
contribution of the fundamental harmonic was averaged out over the whole ensemble of rings 
in the system. The effect of period halving can be easily understood in the free electron 
model assuming that the number of electrons, N , in a given ring rather then chemical 
potential is a conserved quantity. In this case the sign of the rath Fourier harmonics entering 
the expression for the persistent current is (—1)^", and after ensemble averaging (fluctuating 
A^) all the odd terms are averaged out.lll 

The Aharonov-Bohm oscillations of the magnetization with fundamental period $0 were 
first observed in a single gold ring.! The amplitude of the oscillations appeared to be unex- 
pectedly large (on the order of Jq ~ evpjV) and this has triggered an avalanche of theoretical 
papers discussing the AB oscillations in isolated "dirty" normal metal-rings. The key point 
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here is that by properly taking into account the elastic scattering by impurities, the dy- 
namics of the charge carriers in a multichannel ring becomes diffusive (inelastic processes 
are negligible at sufficiently low temperatures). Therefore, a rather straightforward zero 
temperature estimate suggests that Id ~ loiu/L, where in is the electron mean-free path. 
However, the actual value of the current amplitude is greater than by several orders of 
magnitude in the diffusive regime where £d L, and this glaring disagreement has initiated 
a discussion in the literature.iHii 

Metallic rings are intrinsically disordered, generally contain many transverse channels 
and are in the diffusive regime. In this paper we investigate the opposite limit of the effect 
of a single impurity on a single-channel semiconductor ring in the ballistic limit. This limit 
is instructive because it allows a reliable treatment of the particle interactions. Furthermore, 
recent progress in semiconductor quantum wire technology has allowed the production of 
clean low-channel number rings and measurement of their persistent currents.0 Another 
important reason for studying this problem is its connection to the problem of how the 
nature of the Coulomb blockade changes as the number of transverse channels in a Luttinger 
liquid is varied.0 

The measurements of the magnetic moment in impurity-free rings in Ref. revealed 
oscillations with a fundamental period $o and amplitude, Iq. These results are in perfect 
agreement with predictions based on the free Fermi-gas model. At first glance, this result is 
mysterious since at low charge density the Debye screening is weak and long range Coulomb 
forces must be a decisive factor for the quantum dynamics of the system. The effect of 
electron-electron correlations on the persistent current of an isolated metallic ring were 



studied in Ref. 19 where it was shown that in the diffusive regime the interaction can 



enhance the current. For the case of ballistic rings, such a problem was addressed in Ref. 20 



building from the Luttinger model, and in Ref. ^ where the QED radiative corrections 
to the current were calculated. The period and zero-temperature amplitude of oscillations 
are not affected by the Coulomb interactions. This is as it should be according to general 



arguments.cHfE2l Even the temperature dependence of the amplitude is in fact the same as 



for the ideal gas. For an isolated ring with a fixed number of electrons the infiuence of 
temperature is to destroy the oscillations. The exponential suppression of the amplitude of 
the persistent current starts at a crossover temperature, Tq ~ hvp/L, which is of the order 
of the energy-level spacing at the Fermi energy for the free Fermi-gas model. 

In this paper, we study the persistent current in a ID ballistic ring in the limit of 
strong Coulomb repulsion so that the (fiuctuating) one- dimensional Wigner crystal (WC) 
representation of the Luttinger liquid is valid. 0^1 Note that the use of the name Wigner 
crystal here does not imply that the system has long-range order. (It does not.) Within 
a long wavelength approximation the system can be viewed as an elastic chain of identical 
charges completely described by three independent quantities: the (electron) mass m, the 
lattice spacing a, and the sound velocity s. The magnitude of the quantum fiuctuations are 
determinedii by a dimensionless parameter a = -nh/msa. In what follows, we take a -C 1, 
which corresponds to a stiff crystal. 

In the absence of crystal pinning (perfect ring), the persistent current induced by the 
magnetic flux is caused by a uniform sliding of the Wigner crystal as a whole. The elementary 
charge transfer corresponds to a translation of the crystal ring by the lattice constant a. Due 
to the fact that the electrons at the different sites of the Wigner lattice are indistinguishable, 
such a translation in a circular geometry is equivalent to carrying a single charge, e, along a 
closed loop. Therefore the period of the corresponding Aharonov-Bohm oscillations is equal 
to $0) exactly as in the free Fermi-gas model. 

Another simple argument tells us why the oscillation amplitude also remains unchanged 
in the crystalline phase. The uniform translation of the crystal by a lattice spacing a results 
in a "diamagnetic" current Iwc ~ efo/L, and the corresponding "drift" velocity vq ~ Ti/ma 
coincides with the Fermi velocity for a gas of free spinless electrons at the same density. 
Clearly, this very circumstance explains why the cross-over temperature is of the same 
order for the two systems. These simple arguments, of course, can not explain the exact 
coincidence of the zero-temperature expressions in question. This fact follows from a general 
theorem. ii 
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The Wigner crystal in a perfect ring does not get pinned and its persistent currrent is 
(at T = 0) identical to the one carried by free electrons. In order to activate internal modes 
of the chain one has to consider pinning effects. In the latter case, the charge transport 
becomes of the tunneling type (macroscopic quantum tunneling of the Wigner crystal ring) 
and the amplitude of the persistent current acquires a dependence upon both the pinning 
barrier height and the intrinsic properties of the crystal. 

The influence of the barrier on the quantum dynamics of the finite Wigner chain depends 
on the strength of the pinning. Defining = %s/L as the characteristic energy scale in a 
chain of a finite length, L, the strong pinning regime occurs for Vp ^ Ts/a. In this regime 
the quantum fluctuations of the displacement field are effectively cut off at the intermediate 
scale £o ~ CLiTs/aVp) ^ L (see Ref. ^) so that the long-wavelength dynamics [k < vr/fo) is 
quasiclassical and the functional integral for the free energy can be evaluated by the method 
of "steepest descent". 

Charge transfer along the ring can be viewed as a depinning of the crystal by means 
of tunneling. This process has two distinct stages: (i) "fast" tunneling of a small ring- 
segment (/q ^ L) through the barrier, (ii) "slow" relaxation of the elastically deformed 
state developed in the near barrier region. It is the contribution of the relaxation stage that 
dominates the tunneling action.0 Hence, the magnitude of the persistent current crucially 
depends on the elastic properties of the chain. We will show that at low temperatures the 
expression for /($) contains an additional small factor {Ts/ aVpY^°' compared to the case of 
a perfect WC-ring. With increasing temperature the effective tunnel barrier decreases and 
the persistent current starts to grow. Two competing mechanisms — thermally activated 
tunneling and the temperature destruction of quantum coherence due to the enhancement 
of destructive interference — lead (for a <^ 1) to the appearance of a sharp maximum in 
the temperature dependence of the oscillation amplitude at T = — O.ST,. 

Although aX T > Ts the destructive effect dominates and the magnitude of the current 
is exponentially suppressed (this is also true for the ideal crystal), it is partly compensated 
at T ^ aVp by a growing pre-exponential factor {T / aVpY^" . This T^l°' law for tunneling of 
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strongly correlated electrons was predicted in Refs. and 

A relatively weak suppression of the oscillations ensures that at temperatures of the 
order of T ~ aVp the magnitude of the current should coincide with its value for an ideal 
crystal. Here the elastic relaxation stage has no effect since the initial length io{T) of the 
deformed segment is as large as the ring circumference L, and the whole transport is due to 
the uniform tunneling, exp{—Vp/Ts). The same exponential factor determines the amplitude 
of the persistent current in a perfect ring at T ~ aVp. 

With a further increase in temperature, the quantum depinning is replaced by a thermally 
activated depinning process, Vp/Tg — >■ Vp/T, and the corresponding expression for /($) 
acquires at T aVp a factor, exp(— V^/T), that grows with temperature in addition to the 
factor exp(— [7r/2][T/aTs]), which is a decreasing function of temperature. 

The above scenario relies on two assumptions: (i) the crystal is stiff, a <^ 1, and (ii) 
the pinning is strong, Vp ^ Tg/a. In the case of weak pinning the picture of macroscopic 
tunneling of the Wigner crystal ring is changed drastically. Firstly, the long- wavelength 
fluctuations of the displacement field are cut off now but at the scale L (in the weak pinning 
regime C.q is larger than L). The requirement of weak fluctuations puts an upper bound 
on both the chain length and the temperature: L -C ae^^", T < Tg/a. Secondly, for low 
temperatures, T <^ T* (aTsV^)^/^, and weak pinning the charge transport through the 
barrier goes entirely via succesive tunneling rotations of the Wigner crystal ring as a whole by 
a distance a. In accordance with this, the zero-temperature amplitude of the current acquires 



(for Vp ^ To) a small tunneling factor exp(— ay 2m(L/a)lp) = exp(— y/27rV^/To) ^ 1. 
Though at finite temperatures the tunneling rate is increased and the competition between 
the thermally activated tunneling and the enhancement of destructive interference leads to a 
maximum of the persistent current at T ~ T*, this maximum is weakly pronounced. Loosely 
speaking, the effect of moderate pinning, aTg -C ^ T,/a, results in a suppression of the 
zero temperature amplitude of the current and in a shift of the crossover temperature, Tq, to 



the barrier dependent value T* ~ ^oTsVp . For very weak pinning, Vp <^ Tq, the persistent 
current of a Wigner crystal ring is indistinguishable from that of non-interecting electrons. 



At high temperatures, T ^ max{aTs, (aT^l^)^/^}, quantum depinning is replaced by a 
thermally activated depinning process, and the conventional picture of the Aharonov-Bohm 
oscillations in an ideal Wigner crystal is recovered, up to a slight correction due to the 
existence of the barrier. 

To summarize, the temperature dependence of the persistent current differs drastically 
from the predictions based on the Fermi-gas model both in the regime of strong and mod- 
erately weak pinning (Vp > To). To be more specific, the temperature dependence exhibits 
a sharp maximum at T O.ST^ in the case of strong pinning, Vp ^ T^/a. For moderately 
weak pinning, Tq ^ Vp ^ Ts/a, the maximum is weakly pronounced, but its position is 
shifted to a potential-dependent value T* ~ (TqVpY^'^ > Tq which can be considered as 
the crossover temperature (instead of Tq) in this regime of pinning. In both cases the zero 
temperature amplitude is strongly suppressed. This conclusion is supported by the basic 
difference in the tunneling dynamics of strongly correlated and free electrons in an infinite 
one-dimensional system.il'S 

For non-interacting electrons, the influence of a regular potential Vp <^ Ep — where 
Ep is the Fermi energy — is neglig ible.i Although in the presence of scatterers, not only 
thermodynamic quantities but also the energy-levels oscillate as a function of flux, these 
oscillations generally come with alternating signs so that only the energy level closest to the 
Fermi level in effect contributes to the current oscillations. The latter is insensitive to the 
action of a "weak" (Vp ^ Ep) potential. The temperature dependence is also unaffected in 
this case.0 Therefore, measurements of the persistent current in ballistic rings may provide 
a deep insight into the quantum properties of strongly correlated electron systems, and 
hence offer a unique possibility of detecting and investigating properties of a Luttinger 
liquid/Wigner crystal in detail. 
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II. PERSISTENT CURRENT IN A PERFECT WIGNER CRYSTAL RING 



Let us consider a mesoscopic ring formed in a laterally confined two-dimensional electron 
gas. For a low density of electrons {nas <^ 1, where is the Bohr radius) weakly screened 
Coulomb interactions can lead to the formation of a ID Wigner crystalBiE We will study 
the response of such a strongly correlated electron system to an enclosed flux, $, of magnetic 
field directed normal to the plane of the ring. 

It is well knownill that in a non-simply connected space the free energy, F, of the charged 
particles acquires a flux-dependent part resulting in the appearence of an induced magnetic 
moment.il This moment can be related to the persistent current in the ring which is given 



by 



For a non-superconducting ring the amplitude of the persistent current is a decreasing func- 
tion of the ring- circumference, L, and fiux-induced effects can be detected only in samples 
of mesoscopic sizes. 



Adopting the model proposed in Ref. ^ we will describe a ID Wigner crystal as an 



elastic chain of spinless electrons. The input parameters of the model are the crystal period 
a, the electron mass, m, and the sound velocity s. In the continuum limit the Lagrangian 
of such a system in the presence of a pinning potential isS 

^ = ^ - s'i^r} - - cos^)- (2) 

Here (p = 2'Ku{x,t)/a is the dimensionless dynamical field of the crystal {u{x,t) being the 
local displacement of the crystal at the point x and time t) . The notation and ip' indicate 
a time- and spatial derivative, respectively. The pinning potential (with amplitude Vp) is 
assumed to be smooth on the scale of a but well localized on the scale of the ring size L. 
Without loss of generality we place the pinning potential at point x = 0. 

In low dimensional systems one should take into account quantum fluctuations which can 
be very strong in the long- wavelength limit. These fluctuations manifest themselves in the 
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logarithmic increase of the correlator of the displacement field at long distances {x ^ a) and 
lead to the breakdown of crystalline order in infinite Wigner chains even at zero temperature. 
The magnitude of quantum fluctuations in a WC can be characterized^ by a dimensionless 
parameter a 

TT^ /ON 

a = . (3) 

msa 

In what follows we will study a stiff (weakly- fluctuating) WC, which implies that a ^ 1. 

In a perfect (or weakly pinned) WC, long- wavelength fluctuations are cut off at a wave- 
length on the order of the crystal size L. It is physically evident that we can asssume an 
ordered crystalline structure as long as the mean square fluctuations of the dimensionless 
field ip, 

(^')~«/ -IT^oth -—k], (4) 



t/l k ylksT 

are small, so that {ip^) ^ 1 (T is the temperature). For T ^ this inequality imposes a 
restriction in the form of an upper bound on the chain length L <^ ae^/". One can easily 
check that for such samples the thermal fluctuations are suppressed up to the temperature 
T ^ Tg/a {Ts = hs/L). The situation is changed drastically for a strongly pinned WC 
where an "intermediate" cut off scale appears.il The pinned WC will be studied in the next 
section and here we calculate persistent current in a perfect ( Vp = ) Wigner crystal ring. 

In a magnetic field, directed normal to the plane of the ring, the one-dimensional La- 
grangian (^ acquires an additional term, Cint, which describes the Aharonov- Bohm inter- 
actions of the WC with the vector potential of an electromagnetic field, = <I>/L ($ is the 
magetic fiux through the ring). The AB interaction rewritten in terms of the displacement 
field if, takes the form of a total time derivative 

h $ 

and affects, as must be the case, only the quantum dynamics of the crystal. 

The fiux-induced persistent current, Eq. (|I|), is defined in terms of the sensitivity of the 
free energy of the ring to magnetic fiux. For the following analyses, it is convenient to 
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express the free energy, F,, as a path integral over quantum and thermal fluctuations of the 
displacement field. 



ln=-oo •' J 



(6) 



Here Se is the Euclidean (imaginary time) action for the Lagrangian (H), (^, while is the 
total number of electrons in the chain. Since the arbitrary states {f{x, t)} of the WC which 
differ from each other by a constant 27m {n is an integer) are physically indistinguishable, 
one should impose "twisted" boundary conditions on the field ip in imaginary time (see e.g. 
Ref . 0) 

iPnir + P,X) = iPn{T,x) +2-7171. (7) 

This is why the path integration in Eq. (^ includes an additional summation over the ho- 
motopy index n (winding number), which classifies homotopically inequivalent trajectories. 

The appearance in Eq. (^) of a sign-altering weight-factor, which depends on the parity of 
the total number, A^, of electrons in the ring, describes the "parity effect"0i§ in mesoscopic 
systems of strongly correlated electrons. This factor has a simple physical interpretation. 
One can view the important ring-exchange transformations of the WC as a homogeneous 
shift of length a {A(f = 27r) supplemented by (A^ — 1) successive permutations of pairs of 
electrons. The corresponding extra phase 7r(A^ — 1) that appears in the many-particle wave 
function generates the sign-alternating factor in Eq. (P). This factor originates from the 
the fact that the electrons in the WC-ring obey Fermi-statistics and is identical for free- 
and strongly correlated electron systems.ilS This factor does not appear in the infinite 
Luttinger liquid problem@0 but is essential to the physics of finite rings. 

For a perfect Wigner crystal (Vp = 0) the problem of calculating the fiux-dependent part 
of the free energy, AF{ip), can be solved exactly. The Lagrangian is "quadratic" and the 
extremal trajectories obeying the boundary conditions (^ are ^-independent linear functions 
of imaginary time r, 

27m ,^ 
<^n(r,x) = ^r, (8) 
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(/3 = l/Zc^T). By substituting Eq. (||) into Eqs. (|]), (^), and (H) and then performing the 
summation over n, it is easy to find that 

AF{^) = -kBT\nMd,Q) (9) 
g^expflJ), 9 = + 5n (10) 



where d^{9,q) is the Jacobi theta function, Tq = hvp/L (in terms of WC parameters Tq 
aTg, Tg = hs/L), and the parity dependent term 6iy is 1/2 (0) for odd (even). 



By making use of the asymptotic expressions for the "i^a-function (see e.g. Ref. 34) we 
get the desired formulae for the persistent current at high- and low temperatures 



-l)^2/o Jexp (-1 J) sin (27r|^) , T > T( 



(11) 

Jo (1-2 {{1^ + 5^}}), T«To. 



Here {{x}} denotes the fractional part of x. Thus the persistent current carried by an ideal 
Wigner crystal ring is a periodic function of fiux with period $o = hc/e and amplitude 
Jo = evp/L at low temperatures. The current has a paramagnetic character when there is 
an even number of elecrons in the ring (i.e. the induced magnetic moment is parallel to the 
external magnetic field) and diamagnetic for an odd number of electrons. All these properties 
of the persistent current exactly coincide with those calculated using the model of an ideal 
Fermi gas.i The response of a perfect ring (without impurity scattering) to magnetic fiux 
at T = is essentially independent of the character and strength of Coulomb correlations 
among the electrons. For the Luttinger model this was shown first in Ref. |2^. A rigorous 
microscopic derivation of this fact has been given by Miiller-Groeling et al.ii 

We call attention to the numerical factor in the exponent of Eq. ([TlD . According to this 
the temperature, = (2/7c)Tq = (2/7T){hvF/L), that determines the cross-over to the regime 
where the persistent current is suppressed by loss of quantum coherence, I oc exp(— T/Tc), 
is twice as large as for an ideal ring with Fermi-distributed electrondlfl (fixed chemical 
potential fi). The appearence of an additional factor of 2 in the cross-over temperature in 
our case is connected with the statistical propertieSH (thermal averaging over a canonical 
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ensemble with a fixed number of particles) and is weakly dependent on the electron-electron 
correlations. 

Now we need to consider the influence of fluctuations (temperature induced "phonons" 
or plasmons in the WC) on the persistent current of a Wigner crystal ring. If we have a 
rigid displacement of a ring by L/a = N steps (A^ is the total number of electrons in the 
ring) the system returns exactly to its original state, and so do the distortions. For n ^ N 
we must move the distortion to get the system back to its initial state ip{x) (f{x — na). 
This procedure implicitily takes into account the finite period a of the Wigner lattice and 
therefore goes beyond the long- wavelength approximation used in Eqs. (0),@,and where 
all short wavelength processes have been included in the sense that they have renormalized 
the bare parameters of the effective Lagrangian (0). Nevertheless, it seems interesting to 
estimate to what extent our temperature corrections. Eg. (pH]) , are valid. 

To do this it is convenient to go over to momentum space, where the shift (f{x) — <f{x + 
na) corresponds to a simple change of phase in the Fourier transforms of the displacement 
field, (fg — > exp{i9g)ipg, 9g = qna. The action in momentum space reads 



ma 



/ dT{\iPg\'+sV\iPg\'}, (12) 

J —h0/2 



Sll'^h J-hl3/2 

and describes two uncoupled harmonic oscillators with frequency ujg = sq. The extremal 
path corresponds to the classical motion in an inverted potential subject to the boundary 
condition Lpg{T = —hj3/2) = exp{i9g){pg{T = Ti(3/2) . A straightforward calculation of the 
partition function for the model action (|12D gives 



\ , sm\eg/2) 



2M /o^ H-i 

(13) 



smh\(3hujg/2_ 

If we introduce the thermal coherence length = fis/T, then according to Eq. ( pTSj) only 
modes with qC,T ~ 1 are significant and we can use the small-g approximation 

^^e-(We.)^, na<eT. (14) 

Since the number of modes for which qC,T ^ 1 is of the order of L/C,t, we readily get the 
following rough estimate of the action for moving the distortion field by n steps 
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where Ep is the Fermi energy. Comparing Eq.(|T^) with Eqs. (pD-(pA|), it is easy to see that 
the contribution of phonons to the persistent current of a thermally distorted Wigner crystal 
ring is additionally suppressed by a factor a^{T/EpY <^ 1 and can be safely neglected in 
what follows. (It is an irrelevant perturbation.) 

III. MACROSCOPIC QUANTUM TUNNELING OF A WIGNER CRYSTAL-RING 

Now let us study the persistent current in a Wigner crystal ring pinned by a potential 
barrier. A uniform sliding motion of the crystal is impossible in this case and the charge 
transport is inevitably associated with macroscopic quantum tunneling. 

The tunneling of strongly correlated electrons through a weak link in an infinite one- di- 
mensional chain has been studied@0 in the framework of the Luttinger model. It was shown 
that in the case of repulsive interactions the electrons at zero temperature are completely 
reflected by the barrier even for an arbitrarily small height of the barrier. The authors 
of Ref. ^ in essence, came to the same conclusion when studying the tunneling of a WC 
through a potential smooth on the scale a. One of the most surprising results of this paper is 
the claim that the pinning of the Wigner crystal suppresses the long-wavelength fluctuations 
of the displacement field. Unlike the case of a perfect WC, the fluctuations in the pinned 
crystal chain get cut off at the intermediate scale a -C -C which can be found by 
minimizing the energy.!! One finds 



where Vq <^ ms^ is the bare amplitude of the pinning potential. Therefore, even for an 
infinite chain the potential — regardless of its initial strength — is screened only partiallyil 





(16) 





(17) 
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Since the tunneling action of the WC diverges logarithmically at long distances, the charge 
transport along the imperfect infinite ID chain vanishes at T = 0. For a crystal of finite 
size the probability of tunnelling is finite, thus providing a non-zero persistent current in a 
pinned WC-ring subject to a magnetic field. 

In the following we will consider Eq. (|^) to be an effective Lagrangian. According to the 
results of Ref. short-wavelength fluctuations {k < kc) lead only to a multiplicative renor- 
malization [Eq. ([17|) ] of the amplitude of the pinning potential. Therefore the renormalized 
Lagrangian can serve as a starting point for studying the long-wavelength quantum dy- 
namics of a WC-ring. In the limit a <^ 1 (stiff Wigner crystal) one can use quasiclassical 
approximations in calculating the path intergral, Eq. (^), for the free energy. 

Tunneling of strongly correlated electrons — described by the model in Eq. (0) — through 



a smooth barrier has been considered thoroughly in Refs. |20| and In the case of strong 
pinning aVp ^ Tg the process of macroscopic quantum tunneling can be divided formally 
into two different stages: (i) a rapid stage resulting in the formation of a local deformation 
in a small region 2£o (see below) adjacent to the barrier, and (ii) a slow process of elastic 
relaxation of the deformed segment. 

We shall first calculate the tunnel action and persistent current at zero temperature. At 
the rapid stage of the tunneling the displacement field, (p, "hops" {Aip = 2tt) locally in the 



vicinity of the impurity. It is described by the trial instanton-trajectory proposed in Ref. |26. 
The corresponding tunnel action, A^, is of the form 

A. = * (C, + C,^) . (18) 

The exact values of the numerical constants Ci ~ C2 ~ 1 are not essential in what follows; 
2^0 is the initial length of the deformed segment of the WC. It can be found by minimizing 
the total action (see below). 

The second stage of tunneling describes the elastic relaxation of the deformed state and 
the corresponding instanton solution obeys the free equation of motion (in imaginary time) 
for the displacement field ip, 
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Cp + sV" = (19) 

The desired solution should satisfy the boundary conditions 

(^±(r = — oo, x) — (y3±(r = oo, x) = ±27r (20) 
(^±(r,L/2) = <^±(r,-L/2) (21) 

(the subscript, ±, labels trajectories with winding numbers n = ±1) and represent the 
tunneling dynamics for |x| > £o and |r| < tq = £o/s; at the boundaries it has to match 
the trial function describing the first stage of the tunneling. A tunneling trajectory which 



satisfies all these requirements has the form (Fig. la; see also Ref. ^ 



(/9± = TT ± 2 arctan \x\ ' ) '^^^^ 

Here Tg is an arbitrary parameter which denotes the location in time of the instanton. The 
action for the relaxation stage (using the solution (p2|)) is 



It is easy to find that the minimum of the total action At = Ah + Aj. corresponds to the 
following initial length, 2£o, of the deformed segment 

For a stiff crystal this length coincides with the lower bound (|T6|) of the long wavelength 
region, where the effective Lagrangian @) is assumed to be valid. Thus our calculation can 
be quantitatively reliable if Ar':^ Ah. This is the case of strong pinningS aVp ^ Tg, and 

A = ^L(^)+CA. (25) 



a { \ Ts 

The numerical value of the constant C3 in Eq. (^) depends on the first stage of the tunneling, 
which develops in the short-wavelength region, and thus is beyond the accuracy of our 
calculation. 
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With the help of Eq. (^) it is easy to find the fiux-dependent part of the ground state 
energy Eq{^). For a stiff crystal { a <^ 1 ) the tunnel action is large, Af/h ^ 1. Therefore 
we can use in our calculations the well-known instanton dilute gas approximation, where the 
tunneling trajectories are "multi-step" solutions of Eq. ( [T9| ) (Fig. lb; see e.g. Ref. |3^ as well 
as the similar calculation in Ref. ^ of the Aharonov-Bohm effect in charge-density-wave 
conductors) 



Eo = -k^T)im{ Y: E 



n++n- 



■X 



[n±:'oo„+7^=0 ^+'^-' 

g-(ylt/a){n++n_)+j27r0(n+-n_)^^^^_^ |_ ^26) 

Here is the number instantons (anti-instantons), i^' is a numerical constant. At is the 

tunnel action (^, and 6 = $/$o + 6n (see Eq. (|10|) ). 

By doing the sums in Eq. (^) we readily get the well-known general expression for the 
energy of a G-vacuum in the dilute gas approximation (see e.g. chapter 10 in Ref. ^6|) 



A^o('^') = KTs^jAt/he-^'/^ cos 9. (27) 

In the case of strong pinning we only need to keep the main (first) term in the tunnel 
action (^Sj) . Neglecting the inessential prefactors one can derive from Eq. (|27|) the desired 
expression for the oscillating part of the free energy at T = 

AFwc{<^, T = 0) ~ i-l fT, (^-^^ ' e-^^'/" cos (27r^) . (28) 

We see that the amplitude of the persistent current in a pinned WC-ring is much smaller 
than in a perfect crystal. It decreases rapidly with increasing crystal size ~ (1///)^+^/°". Such 
a dependence for a <^ 1 is intermediate between the regimes of Aharonov-Bohm oscillations 
in metalsi (/ ~ 1/L) and insulatorsi'i (J ~ exp(— L/£g), ig = hs/Ag, Ag being the gap in 
the energy spectrum at the Fermi level). 
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IV. ANOMALOUS TEMPERATURE DEPENDENCE OF A PERSISTENT 
CURRENT IN A STRONGLY PINNED WIGNER CRYSTAL 



In a perfect Wigner crystal the temperature always acts to destroy the Aharonov-Bohm 
oscillations (as in the case of a free Fermi-gas). Starting from the cross-over temperature 
To ~ hvp/L the oscillation amplitude is exponentially suppressed. In the presence of a pin- 
ning potential the charge transport is due to tunneling and the influence of temperature can 
be more complicated. From the argument that a finite temperature facilitates tunneling, one 
can expect an increase of the persistent current in the low temperature regime. Combining 
this effect with the loss of phase coherence at high temperatures due to the enhancement 
of destructive interference, we may spectulate that the temperature dependence of the per- 
sistent current exhibits a maximum at T ~ = /is/L. Below we will show that for a stiff 
crystal this maximum develops into a sharp peak. 

To prove this statement we have to calculate the flux-dependent part of the free energy 
of a WC-ring at non-zero temperatures. We will again divide the process of depinning by 
tunneling into two stages: (i) a rapid formation of a distorted region of size -C L in 
the vicinity of the impurity, and (ii) a slow relaxation of the elastic tension during a "time" 
Tp = hp {f3=l/T). 

In the previous section it was shown that in the case of strong pinning, aVp ^ Tg, the 
dominating stage is the process of elastic relaxation of the deformed state. Therefore we 
will not take into account the temperature corrections to the action Ah, leaving Eq. (p!8D 
intact (now, however, the initial length of the relaxing segment 2ij3{T) becomes temperature 
dependent). 

The main temperature dependence of the total tunneling action comes from the relax- 
ation process due to the deformation of the Larkin-Lee instanton, Eq. (0). As before, we 
seek the solution of the free equation of motion, Eq. ([T9|), but with twisted boundary con- 
ditions imposed on the flnite interval of imaginary time [0,/5] (see Eq. (|^) for n = ±1). The 
desired solution should transform into the Larkin-Lee instanton at low temperatures and 
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describe a uniform sliding of the WC-ring, Eq. in the high temperature region. 
It is easy to verify that the exact solution is given by 



(faix, t) = n + 2a arctan 



, /7rb|\ f 71 f hp^ 

coth - — - tan — — r 

\hsf3l \hf3\ 2, 



(29) 



where a = ±1. This "periodic" instanton solution holds when the spatial coordinate, x, 
is outside the interval [—^13,^/3] containing the pinning potential, and for imaginary time 
\t — hp/2\ > i/^/ s. It should be noted that our solution, Eq. (^) is similar to the so called 
"colorons" in quantum chromodynamics0, and it can be found using a method proposed in 



Ref. 41. 



In the low- and high-temperature regions Eq. (^) takes the following asymptotic expres- 
sions 

~ TT ± 2 arctan (^ ^^ ~\Jf^^ ^ ' ^ °° (^0) 

Svrr 27rr ^ , , 

, ^-^2.-—, P^O. (31) 

(at low temperatures the interval of imaginary time is infinite, j3 00, and the center 
of the instanton can be placed at the arbitrary point r^). Equation ( pO]) is nothing but 
the Larkin-Lee instanton.ll The high temperature limit, Eq. (pT]), describes a homogeneous 
instanton (v?+) and anti-instanton (v9„) for the winding numbers n = ±1. Notice that 
at finite temperatures the center of the periodic instanton, Eq. (P9|), is a fixed quantity 
rather than an arbitrary parameter (it is placed exactly in the middle of the imaginary time 
interval) . 

The action for the tunneling trajectory ( ^91) is 



Ap = - lln 

a 



sinh(7rT/Ts 
sinh(27rr/T/ 



(32) 



where Th = hs/i/3{T). The initial length, 2£^(T), of the deformed segment which minimizes 
the total action is equal to 

^ 27rT V2 27r T / ^ ^ 
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We will assume that ifs L, a criterion which one can show to be equivalent to a restriction 
on temperature, T <^Vp. 

By substituting Eq. (^) into Eqs. (|18|) and ( |32D one finds for the total tunneling action 



--— + lnsmh (^TT— ) + 2 In 



^ ^ 27r T . 
ChaVp 3 
27r T 2 



27r T 



(34) 



Here Ci ~ C2 ~ 1 (see Eq. (0)) and Eq. ( |3^ is valid at temperatures T -C aVp. The 
tunneling action, considered as a function of temperature, has a minimum at r — ~ 
0.5Ts, which is independent of the values of the numerical constants Ci^2 {Tm is a solution of 
the transcendental equation coth(2;) = 1/2 + 1/z, where z = ttT/Ts). Thus it characterizes 
the relaxation stage of the macroscopic quantum tunneling of a Wigner crystal ring. The 
extremal action, At(Tm), divided by the tunnel action at T = 0, Eq. (p5D , is also entirely 
determined by the relaxation stage. The minimum of the action corresponds to a maximum 
in the persistent current and hence the temperature dependence of the AB-oscillations in a 
strongly pinned WC-ring turns out to be anomalous. 

In the quasiclassical approximation (stiff WC: a <^ 1, strong pinning: aVp ^ Tg) the 
trajectories with smallest winding numbers (n = ±1) give the main contributions to the 
oscillating part of the free energy. The persistent current in this case is determined by its 
fundamental harmonic 



Iwci^.T) = i-lflwc{T)smi27r $/<l>o) . 



(35) 



The amplitude of the current calculated for the "periodic" instanton solution, Eq. (|^), is 



T 

Iwc{T) = lwc{0)—exp 

J- .'f 



'-fa 



T 

1 s 



aVr, 



l/a 



<T<^aVp 



71 



f(x) = -X - In 



sinh(7rx) 



nx 



(36) 



where Iwci^) is the value at which the persistent current saturates at zero temperatures 
(due to the finite length of the ring providing an infrared cutoff scale T^): 
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l/o 



According to Eq. ( pGf ) the persistent current for a ^ 1 peaks (see Fig. 2) at T = T, 



O.ST^. The width of the peak is of the order ^/T^Tg. Even though the oscillation amplitude 
at the maximum is much less than the value of the persistent current in a perfect ring 
at the same temperature T^, the temperature- induced suppression of the AB-oscillation is 
comparatively weak since in the regime T ^ Tg 

This is because of the thermal stimulation of tunneling, which gives rise to the monotonically 
increasing prefactor in Eq. (P^). The unusual temperature dependence — a power-law 
behaviour, T^/" — of transport properties of strongly correlated electrons was discussed 
previously in Refs. ^ and 

Rigorously speaking, the expressions derived in this section are valid at T <^ aVp. Nev- 
ertheless, we can use them for a qualitative analysis of the situation at T ~ aVp. According 
to Eq. (|33|) the initial length of the deformed segment in this case is of the order of the 
ring circumference (2£^ — > L at T ^ C2al^/vr) and the relaxation stage of the tunneling 
vanishes. Now a shift of the WC-ring by tunneling a distance a corresponds to a uniform 
motion of the crystal as a whole. We can therefore expect that the persistent current will 
coincide with that of a perfect ring. In fact, according to Eq. (|3^), at T ~ aVp 

Iwc{T ~ aVp) ~ exp (^-const.^^ . (39) 



This expression coincides (up to a numerical constant) with the corresponding result (|TT]) for 
the perfect WC-ring. With a further increase of temperature, the pinning can be removed 
by thermal fluctuations and the AB-oscillations take place in a weak pinning regime. 



V. PERSISTENT CURRENT IN THE WEAK PINNING REGIME 



So far we have studied the coherent properties of a Wigner crystal ring assuming that 
— despite the "short-wavelength" {k < kc) screening of the bare potential Vq, — the renor- 
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malized height of the potential barrier Vp, Eq. ([T7|), is large enough to classify the tunneling 
dynamics of the WC-ring as being in the strong pinning regime, where Vp ^ T^/ a ~ ms^a/ L. 
In this case the coherent properties of the WC are radically different from those in a perfect 
ring up to temperatures T ^ aVp. Now we will concentrate on the influence of weak pinning 
{Vp -C Tg/a) on the magnitude of persistent current. 

In the case of weak pinning, the influence of the barrier on the intensity of quantum- 
and thermal fluctuations of the displacement field in a WC of a finite length, L, is negligibly 
small. Because of this we can invoke the same restrictions on the allowed length of the rigid 
chain and the temperature interval as in the case of the ideal crystal: L <^ ae^^". With 
these limitations, we can neglect the renormalization of the height of the potential barrier, 
Vr = V^exp(— i < (/}^ >) ~ Vp{L/a)°' ~ Vp, and the Lagrangian (||) can be regarded as a 
quasiclassical model of a stiff (a <^ 1) Wigner crystal. 

In the regime of weak pinning it is reasonable to use perturbation theory for calculating 
the persistent current. The equation of motion in imaginary time for the dynamical field, 
(fi, takes the form 



In zeroth order the tunneling trajectories are coordinate independent functions (pn{T), which 
are linear in imaginary time, Eq. (§). Hence, the first order correction to the Euclidean action 
of an ideal WC is 



This action describes the ordinary thermally activated charge transport "over the barrier" . 

In order to find the next order corrections we will look for a perturbation solution of 
Eq. ( P0| ) of the form r) = C(j„r + r), where Un = 2'Kn/h(3 and the function 

ipn{x,T) satisfies periodic boundary conditions 



(fi + s^ip" — Vpb{x) sin ^9 = 0. 



(40) 



A^'V^ = VpiT 



(41) 



V'n(a;, r + /3) = -^Jyx, r), ?/'„(0, r) = ^/'„(L, r). 



(42) 



It is easy to find that in second order i 



IS 
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(2), ^ Vp sinh(a;„x/s) - sinh(a;„(x - L)/s) . 

i^ni^^-r) = TT — rj-. : smKr) 

ZuJnS cosn(tt;„L/s) — 1 

. [ K \^ sinh(u;„L/g) [sinh(2u;„L/g) - sinh(2u;„(x - L)/s)] . 
Uujns) (cosh(cu„L/s) - 1) (cosh(2cu„L/s) - 1) ^'""^ ^ ^ 

By substituting r) into the Lagrangian (|^) and performing the spatial and (imaginary) 
time integrals, one finds the Euclidean action, S'^, of a weakly pinned WC to be 

r, ,^ TT T 2 ^ K 3acoth(7rnT/T,) ,,,, 
Se n = + n2m — + ^ + — ^ '—^ (44) 

(note that in the last term of Eq. (^4]) n 7^ 0). 

The desired expression for the oscillating part of the free energy is found by substituting 
the action (^) into Eq. (^) and performing the summation over the homotopy index n. It 
is easy to see from Eq. (^) that at T — > the contributions of the terms stemming from 
the pinning potential exceed the action of the ideal ring and thus the perturbation theory 
fails. In the low temperature region T <^ (aT^V^)^/^ the calculation of the persistent current 
again requires non-perturbative methods and will be carried out below. 

We consider the most interesting case of moderately weak pinning, when the barrier 
is still high enough that aTg ^ Vp <^ T^/a. (For still weaker pinning the corrections to 
Eq. (11) are negligibly small.) It is easy to see that for (moderately) weak pinning the 
tunneling of the WC-ring should be homogeneous. This is because any distortion of the 
crystal structure results in an effective barrier height Vei ^ ms'^{a/ L) ~ Tg/a, which by far 
exceeds the applied potential Vp. Thus the trajectories which describe tunneling in the weak 
pinning case are independent of spatial coordinate and as a function of imaginary time have 
the form illustrated in Fig. lb. They consist of elementary tunnel rotations (A(y9 = ±27r), 
occurring at arbitrary points in imaginary time and in arbitrary sequences subject to the 
only restriction that the net rotation in the "time" interval [0, (3\ equals 2tt. 

A single- instanton (anti-instanton) trial trajectory, V'+(-), of the desired form 

(^+ _ = ±2tx— , (45) 

where the parameter <^h/T defines the time scale of the tunneling rotation (see below), 
can be justified as follows. It is easy to verify that the one-parameter trajectory 
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7~ (yV Stt t 

ip{T,x;To) =271 |-To^-^sinh( — |x|)sin(27r — ), (46) 

To h TqS To 

where tq is an arbitrary parameter, is an exact solution of the nonhnear Eq. pO|). The action 
for this set of trajectories takes the form 

A{to) = ^ + ToVp + —^T^ smh 27r— . 47 
2ro oTl tqs 



For weak pinning the last term in Eq. p7|) can be treated as a perturbation and the 
extremal action is easily found by an iteration procedure. After some strightforward algebra 
we get 



'i^; (48) 



At{n) = 27r^ = a^2m^Vn, (49) 

where the effective barrier height, Vr, is 

Vn = V,il + l^). (50) 

For weak pinning, aVp -C Tg, the renormalization of the bare potential, Vp, is small thus jus- 
tifying the use of the trial instanton trajectory ( ^5] ) and the use of a perturbation expansion 
in solving the equation for the minimum action. It is worth noting that the last expression 
in Eq. ( ^9]) shows explicitly that the process studied describes a homogeneous tunneling of 
a Wigner crystal ring over a distance a. 

Since in the low temperature region, where T <^ T* = (aVpTsY^"^ , the characteristic 
time, T^, of an elementary tunneling rotation, Aip = ±27r, satisfies the inequality Tr^, <^ h, 
one can regard the successive rotations as uncorrelated events and may use an instanton 
dilute gas approximation when deriving the fiux-dependent part of the free energy. Thus 
the instanton described by Eqs. (|i5| ) and (|i8| ) plays the same role in the weak pinning regime 
as the Larkin-Lee instanton, Eqs. ( p2D and (^), does in the strong pinning regime. From 
Eqs. (23) and (24) we immediately get the zero-temperature limit of the persistent current 
{aTs ^ V^p < Ts/a) as 
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(51) 



Notice the exponential supression of the Wigner crystal current by the impurity potential. 
This is in sharp contrast to the case of noninteracting electrons, where even a strong potential 
(of the order of the Fermi energy) only leads to a power-law supression.@ 

The qualitative picture described above enables one to speculate about the temperature 
dependence of the persistent current in the weak pinning regime. Since temperature helps 
tunneling, we may expect an anomalous temperature dependence of the persistent current 
of a weakly pinned Wigner crystal ring. 

It is physically evident that at finite temperatures the effective barrier for tunneling is 
smaller than at T = 0. One can easily take this effect into account by making use of a 
one-particle analogy of Wigner crystal tunneling. In our case (see Eq. (^)) the stiff crystal 
with a total mass mN will tunnel through the barrier at T 7^ via thermally excited states. 
The corresponding temperature dependent action. At, for T <^ T* can be estimated as 



follows (see e.g. Ref. |25|) 



-Ar/h 



(here we neglect the distinction between Vr and Vp). 

At low temperatures the upper limit of integration in Eq. (|52D is irrelevant and we readily 
get the finite temperature corrections to Eq. (HO) as 



T < Vp. (53) 

The temperature corrections to the tunneling action — which are linear in T at low temper- 
atures, much as in the case of strong pinning, Eq. (30) — diminish the action and thus lead 
to an enhancement of the persistent current. The increase of persistent current amplitude 



saturates at a temperature T ~ T* = y^VpTo when the competing temperature- induced 
loss of quantum coherence becomes important. In contrast to the case of strong pinning the 
difference between the tunneling actions at the mimimum and at zero temperature is not 
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large {^At/h ~ 0(1)). The maximum in the amphtude of the persistent current is therefore 
weakly pronounced. Thus the only distinctive feature of the temperature dependence of a 
moderately pinned {Vp ^ To) Wigner crystal ring — as compared to the ideal Fermi gas — 
is a shift of the crossover temperature Tq ^ T* ^ Tq. For very weak pinning, Vp -C Tg, 
even this distinction vanishes and the response of a Wigner crystal ring to magnetic flux 
approaches that of non-interacting electrons. 



The high-temperature region, T > maxjaT^, JaTsVp}, can be described by perturbation 



theory, Eq. (|4^) . In this case the main contribution to the free energy comes from the 
trajectories with smallest winding numbers, n = ±1. Hence, according to Eqs. ( ^If ) and @, 
one finds 

A.,*) . 2r,-l)~e.p {^l^y^-'i coth (4) (^)] 00. (2.^) , (64, 
We see that at high temperatures weak pinning leads to additional suppression of the cur- 



rent in a perfect ring, Eq. ([Tl|), and this potential-dependent contribution decreases with 
increasing temperature. 
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FIGURES 

FIG. 1. Sketch of tunneling trajectories of the dynamical displacement field ipn{T) of the Wigner 
crystal-ring. The trajectories in (a) correspond to 'single-step' tunneling and winding number n = 1 
(full curve) and n = — 1 (dashed curve). In (b) an example of a n = 1 'multiple-step' tunneling 
trajectory is given. 

FIG. 2. Temperature dependence of the normalized persistent current in a strongly pinned 
Wigner crystal of different stiffness (measured by = 2msa/h; Tg = hs/ksL, sec text). The 
sharp peak for stiff crystals is a result of a competition between two effects of temperature: a 
reduced renormalized tunneling barrier and a loss of quantum coherence. 
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